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Abstract 

We consider non-standard matter effect in flavor conversion of neutrinos crossing 
the core of the Earth. We show that oscillation of core crossing neutrinos with 
E > 0.5 GeV can be well described by a first order perturbation theory. We show 
that due to non-standard matter effect varying chemical composition in the Earth 
can modify the neutrino flavor conversion by 100%. Effects of CP violating phases 
in non-standard Neutral Current interactions are emphasized in particular. 

PACS: 14.60.Pq, 13.15. +g 

1 Introduction 

It is well known that non-standard interaction(NSI) can induce non-standard matter effect 
for neutrino oscillation in medium. Neutrino flavor conversion induced by non-standard 
matter effect was proposed as a candidate solution to the solar neutrino anomaly pQ. The 
present experiments told us that LMA MSW solution [H [2] with the standard interaction 
is the solution to the solar neutrino problem [3], IH EJ [6] . Non-standard matter effect is 
small in oscillation of solar neutrinos. However non-standard matter effect can be much 
larger for neutrinos with high energy (E > 10 GeV), e.g. for long baseline neutrinos, 
atmospheric neutrinos, cosmic neutrinos from the galactic or extra-galactic sources. This 
is because flavor conversion induced by flavor mixing in vacuum decreases as energy 
increases while the matter effect does not decrease with energy. Previous works on effect 
of NSI in neutrino oscillation include [7J [S] . 

Non-standard matter effect can be induced by non-standard Neutral Current interac- 
tion of neutrinos with electron, proton and neutron. In this respect non-standard matter 
effect in neutrino oscillation is not only a way to probe physics beyond the Standard Model 
but also a way to probe chemical composition in matter. Incorporating non-standard mat- 
ter effect in neutrino oscillation introduces more CP violating phases in the Hamiltonian. 
These CP violating phases interfere with the CP violating phase in vacuum and can give 
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interesting phenomena. In matter with varying chemical composition these CP violating 
phases can contribute with different combinations in observables. 

It is the purpose of the present article to study the effect of varying chemical compo- 
sition in the Earth. Effects of CP violating phases in the non-standard interaction will be 
analyzed in particular. In section [2] we show that oscillation of core crossing neutrinos in 
the Earth can be well described by a first order perturbation theory which was developed 
in a previous paper by the author. Scenarios with different CP violating phases and vary- 
ing composition in the Earth are shown. In section [3] we show the effect of non-standard 
interactions and CP violating phases in the non-standard interactions. We summarize 
and comment in section HI We do analysis using the density profile of the Preliminary 
Earth Model(PREM) [9]. 



2 Non-standard matter effect in the Earth 

We consider oscillation of three flavors of neutrinos: ip = (v e , v^, v T ). The evolution 
equation is 

i—ip(x) = H(x)v/j(x), (1) 

(jjJb 

where 

H(x)=H + V(x), (2) 

H Q = ±U diag{0, Ami, ^ . (3) 

V(x), a 3 x 3 matrix, is the potential term accounting for the matter effect. U is the 3x3 
neutrino mixing matrix in vacuum. U is parameterized using standard parameters 612, 
#13, #23 and #13, the CP violating phase. 

In the presence of non-standard NC interaction the potential term can be written as 
follows 

/ o v efl V eT \ 

V(x) = diag{K, 0, 0} + V, e V^ T , (4) 

\Vre V T , V TT ) 

where V e = \f2GpN e is the potential with standard charged current interaction. Gp is 
Fermi constant and N e is electron number density. Vm is from non-standard NC inter- 
action. V t l = Vki because the Hamiltonian is hermitian. x dependence in Vki has been 
suppressed in Eq. (J3J). V ee has been made zero in our convention. This is achieved by 
shifting the phases of neutrinos: V\ — > e~ l f dx Vee v\. 
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In this convention Vki is 

V kl = V2G F J2 (/« - fee)Ns 



s=e,p,n 



Ve [ (/« ~ fee) + Ukl ~ f$Rn], , (5) 



s=e,p,n 



where 



R n = (N n - N e )/N e . (6) 



fkii fki anc ^ fkl are ^ ne dimensionless strengths of non-standard four Fermion interactions 
v^2 /jy Gi? s7 M s Vkl^Ui- N p and iV n are number densities of proton and neutron in matter. 
In obtaining the second line of Eq. §5§ N e = N p in neutral matter has been used. 
We can re- write V(x) as 

e^e e Mr . (7) 

fre 6 rr y 

where = V^/K. We can write 

e k i = 4iO- + R«r kl e- i + u ), (8) 



where 



c fcz - 

s=e,p,n 

is constant in matter, depends on the chemical composition in matter and may 
have x dependence in neutrino trajectory. Tki and <fiki are real numbers. e° ee = and 
r ee = in our convention, = and 0/^ = — (fiik because of the hermiticity of V. 

Constraints on come from direct test on NSI [TUl [TT] and the neutrino oscillation 
experiments. Test on NSI can not be directly translated to constraint on e^. These 
constraints have been discussed in our previous work [12]. It was shown that present 
constraints are |e MAt |, |e TT | < 10~ 2 and |e Me |, |e MT | < 10~ 2 , |e er | < 10 _1 [HI ED5 EE] . 

It is clear that e° kl introduces three CP violating phases in addition to the phase 5\z 
in matrix U. They are phases of e° r „ T . Furthermore ^e^er^r become independent 
phases in case that chemical composition varies in matter. So in matter with varying 
chemical composition, i.e. R n not a constant, we have seven physical CP violating phases 
in total. They can give interesting phenomena in neutrino oscillation. In Earth matter 
R n is estimated[15j 

fO. 024, mantle , , 

n ~ 1 0.146, core 1 ' 
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Figure 1: Left Piy^ — ► u e ) versus energy; right P(y e — > u T ) versus energy. L = 12000 km, 
Amlj = 8 - x 10_5 eV 2 , Am§ 2 = 3. x 10" 3 eV 2 . sin 2 2fl 23 = 1, tan 2 12 = 0.41, sin 2 20i 3 = 0.01, 
<fr 13 = tt/6. = 0.01 e" i7r / 20 , e° er = 0.04 e" i7r / 3 , e° T = 0.01 e"^/ 20 . r = 5, <j& e/1 = <^ T = vr/2, 
er = 0. PREM density profile is used for computation in this figure and all remaining figures 
in this article. 

We will use numbers in ffTUl) in our analysis in the present article. We will concentrate on 
neutrinos with core crossing trajectories. 

It was shown in a previous work that oscillation of neutrinos in the Earth can be well 
described by a first order perturbation theory [ 1 2], [T6] . The theory was analyzed with the 
assumption that is a constant in neutrino trajectory. We show in this section that this 
theory works perfectly well taking into account the fact that chemical composition in the 
core and in the mantle are different. 

We quickly review the perturbation theory. We denote L as the length of neutrino 
trajectory in the Earth. For core crossing neutrinos(L > 10690 km) we write the evolution 
matrix M as 

M = M 3 M 2 M U (11) 

where M 2 is the evolution matrix in the core and M 13 are evolution matrices in the mantle. 
< x < Li and L 2 < x < L are the parts of trajectory in the mantle; L\ < x < L 2 is 
the part of trajectory in the core. We average potential in the mantle and in the core 
separately 

Vi= 1 ^ dxV(x), 2 = 1,2,3, (12) 

where L% = L. Using Vi we define the average Hamiltonian 

+ % i = l,2,3 (13) 
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Figure 2: Time reversal asymmetry, At, versus energy. Left r = 5; right r = 2. Other 
parameters are the same as in Fig. [TJ 



and get eigenvector and mixing matrix 

HU m i U m i ^ 1) 2, 3. 

Aj is a vector. The evolution matrix Mj is expressed as 

Mi =U mi e- i & Li - Li -*\l-iC i )Ui li , i = l,2,3 
C l is a 3 x 3 matrix accounting for the non-adiabatic transition: 



where 



dx e^ x Uli 5Vi(x) U mi e-^ x , 



5V i {x) = V{x)-% 



(14) 



(15) 



(16) 



(17) 



It is clear that (C7*)t = C* holds. 

In [12] we have discussed in detail that this theory is indeed doing expansion using 
small quantities. Cj fc (j ^ k) is suppressed by small quantities for neutrinos with E > 0.5 
GeV. Second order effect is of order 0(C 2 ) and is further suppressed. 

In Fig. [T]we compare the result of numerical computation with that computed in the 
first order perturbation theory, i.e. using Eqs. (fTTj) and (j!5p . For simplicity we have set 



0, r eil 
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We see that result computed using the perturbation theory is in remarkable agreement 
with the numerical result. 

We also show the zeroth order result, i.e. result computed by setting Cj zero in Eq. 
ffTTl) . The zeroth order result is an analytical result computed using average potentials 
in the core and in the mantle separately. We see that the analytical result is not a bad 
approximation to the oscillation pattern. It qualitatively describes the neutrino oscillation 
pattern and can help a lot when making qualitative discussions. 

In Fig [2] we show plot of time reversal asymmetry versus energy. At is defined as 

A = P(Ve -> ~ P(V» -> Ve) 

T P( Z , e ^^)+P(^^ Z y e )- 

Again we see that the first order perturbation theory gives a perfect description of the 
oscillation pattern. The analytical result gives a qualitatively good approximation to 
neutrino oscillation. It can help in making qualitative discussions. 



(19) 



3 Flavor conversion of core crossing neutrinos 

In this section we illustrate the effect of CP violating phases of NSI in neutrino oscillation. 

As shown in the last section, oscillation of core-crossing neutrino can be qualitatively 
described by approximation which uses average densities in the mantle and in the core 
separately. This is an analytical description. We use this description to simplify the 
discussion and see the effect of 4>ki in neutrino oscillation. 

It is easier to discuss in the large energy region where we can re-write the Hamiltonian 

as 

H = V + H h (20) 

where 



Vb = diag{V e ,0,0}, H X = V-V Q + H Q . (21) 

Vo is taken as the leading term in the Hamiltonian. Hi is taken as perturbation. V — Vq 
is for the non-standard matter effect and Hq is the Hamiltonian in vacuum. Hq decreases 
as energy increases. 

Using the average potentials in the core and in the mantle we can get the evolution 
matrix using perturbation in H x . As an example, v e — > v T amplitude is 

A(u e - v T ) « { -^(e-^ - 1) + ^(e-^ - l)e*» 

' * €3 



+ /m ( e ~ iV1 ~ l)e~^ c+V5l) , (22) 
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Figure 3: P(y e — > v T ) versus energy, L = 12000 km; P{u e — > v T ) versus distance L, E = 50 
GeV. 4>en — <pfvr — 7r/2. Other parameters are the same as in Fig. [TJ 

where p c = V e c (L 2 — Li) and (pi = V™L\. H™ and H° are averages of Hi in the mantle 
and in the core. V™ and V e c are averages of V e in the mantle and in the core. The property 
of approximately symmetric density profile in the Earth has been used in Eq. (1221) . It 
can be written as 

i fJ m \ fJ c j-f m 

A{v e - Vr ) * ^^(e-- - 1) + (-i - ^U(e~^ - l)e--, (23) 

where if = tp c + 2<pi. 

Neglecting terms of order 0( ^™y , sin 6*13), we get 

A{v e -> v r ) w e° e (1 + 0.024 r re e^ e ) (e - ^ - 1) 

+0.122 e° e r re - l)e~ iipi . (24) 

Eq. (jTUJ) has been used in obtaining Eq. (1241) . A(i/ e — > z/ r ) is determined by e° e modulated 
by contribution of r re and Te . For neutrinos which do not cross the core of the Earth 
the second term in the r.h.s. of Eq. ([24"]) is absent. 

One can see clearly in Eq. (1241) that if r Te = the transition amplitude, A[y e — > i/ T ), 
is determined by e° e and modulated by factor e _lip — 1. Hence P er is proportional to 
function sin 2 (<p/2). In the right panel of Fig. [3] we see plot for this case. For r re = and 
4> Te = 0, P eT has three peaks with roughly equal heights, as expected. For r re = 5, P er is 
considerably changed. When {p re = and ip re = 7r, P eT is considerably modified around 
the third peak (for core-crossing neutrinos). The height in this peak is quite different from 
that in the first peak (for neutrinos crossing the mantle only). This is quite different from 
the case with r re = 0. In the left panel of Fig. [3] we show the plots of transition probability 
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Figure 4: P(u^ — > u e ) versus energy, L = 12000 km; Piy^ — > z^ e ) versus distance L, E = 50 
GeV. (/> e ^ = <^ T = 7r/2. Other parameters are the same as in Fig. Q3 

versus energy. Again we see the effect of r re and y? re . P er is also slightly modified in the 
first peak when r re ^ 0. This is because of the correction by r Te to the Hamiltonian in the 
mantle, as shown in the first term in the r.h.s of Eq. (|24|) . Comparing with the transition 
probability in neutrinos crossing the mantle, the core-crossing neutrino events encode the 
information of r re and (p Te . And effect of ej^ and r^i are distinctly different in oscillation 
probability. 

In Fig. [3] we see that when L = 12000km P eT is reduced when (p re = and is enhanced 
when (p Te = it. This phenomenon can be understood by considering an interesting case 
which happens when 

+ <Pi ~ 2n7r, (25) 

where n is an integral. Hence (p ~ (pi + 2rm. This is the region of parametric resonance 
for oscillation with standard matter effect [TBI \W\ . In the presence with non-standard 
matter effect we see that the amplitude is not always enhanced. Using Eqs. fl24|) and (f2~5l) 

we get 

A[v e -> v T ) w e° e (1 - 0.098 r re e~ i(t,Te ){e~ iLp - 1). (26) 

The amplitude is reduced for <ft Te = and is enhanced for <p Te = ir. When r = 5 the 
transition probability is reduced or enhanced by 100%. When er = 7r/2 P eT is not much 
enhanced. This is understood by noting that according to Eq. ( 1261) P er is enhanced by 
factor 1 + (0.098r) 2 w 1 + O.Olr 2 . It is a 25% increase when r = 5. In Fig. H] we plot 
-P^a* — > u e) versus energy and the distance L. We can also see the effect of and <pki m 
this plot. In the right panel of Fig. H] significant modifications are seen in the second and 
third peaks. 
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4 Conclusions 



In summary we have analyzed non-standard matter effect in ffavor conversion of neutrinos 
crossing the core of the Earth. We have shown that a first order perturbation theory gives 
a perfect description of neutrino oscillation for core-crossing trajectories. The analytical 
description, which uses only zeroth order result, gives a good approximation. 

One interesting thing is that there are six physical CP violating phases associated 
with the non-standard matter effect when chemical composition changes in matter. This 
is what happens to core crossing neutrinos. It is different from the case when chemical 
composition does not change. In the latter case there are only three physical CP violating 
phases. We analyze effect of additional CP violating phases in neutrino oscillation. 

We have shown that due to non-standard interaction different chemical composition 
in the core and the mantle ( different N n /N e ) can modify neutrino flavor conversion by 
100%. We analyze in particular the region of parametric resonance. It is shown that 
in this region the non-standard matter effect does not always give enhancement to the 
amplitude. Depending on the CP violating phases the non-standard matter effect reduce 
or enhance the neutrino flavor conversion. The signature of non-standard interactions lies 
in the dependence of the neutrino flavor conversion rate on E, the energy of neutrinos, 
and L, the length of neutrino trajectory in the Earth. To figure out these interactions we 
need neutrino sources with different energies and baselines. 

The analysis presented in the present article shows that core crossing neutrino events 
provide an interesting way to test interactions of neutrinos beyond the Standard Model. 
They also provide an independent way to test chemical composition in the Earth. 

Acknowledgment: The research is supported in part by National Science Foundation 
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